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Derivation  of  approximate  formulas  for  determining 
the  added  mass  and  damping  coefficients  of  two- 
dimensional,  small-waterplane-area,  twin-hull  (SWATH) 
sections  is  described.  The  added  mass  and  damping  co- 
efficients of  interest  are  those  associated  with  a 
forced  oscillation  of  SWATH  sections  in  heave,  sway,  or 
roll  mode  in  a free  surface.  The  objective  of  deriving 
the  approximate  formulas  for  the  hydrodynamic  coeffi- 
cients is  to  simplify  the  computation  of  motion  of 
SWATH  ships  in  waves  without  sacrificing  its  accuracy 
significantly. 

The  approximate  formulas  are  derived  based  on  the 
potential-flow  theory.  The  damping  coefficients  are 
obtained  in  terms  of  the  outgoing  wave  amplitudes,  and 
the  added  mass  coefficients  are  obtained  by  using  the 
damping  coefficients  through  the  so-called  Kramers- 
Kronig  relations.  Within  the  frequency  range  of  practi- 
cal interest,  the  approximate  formulas  provide  satisfac- 
tory results. 

ADMINISTRATIVE  INFORMATION 

This  work  was  sponsored  by  the  Naval  Material  Command  as  part  of  the 
High  Performance  Vehicle  Hydrodynamic  Program  of  the  Ship  Performance 
Department,  David  W.  Taylor  Naval  Ship  R&D  Center  (DTNSRDC).  Funding  was 
provided  under  the  Ships,  Subs  and  Boats  Program,  Element  62543N,  Task 
Area  ZF-43421001,  Project  Number  ZF-43421,  Work  Unit  1-1500-102. 

INTRODUCTION 

One  of  the  advantages  expected  from  a small-waterplane-area,  twin- 

hull  (SWATH)  configuration  is  the  improved  seakeeping  qualities  compared 

to  monohull  ships  in  moderate  sea  conditions.  To  assist  in  the  development 

of  the  SWATH  concept,  an  analytical  prediction  of  motion  of  SWATH  ships  in 

waves  has  been  developed  at  the  David  W.  Taylor  Naval  Ship  R&D  Center 
1*  2 

(DTNSRDC).  ’ The  analytical  method  can  predict  the  motion  of  a SWATH 
ship  in  five  degrees-of-freedom  in  regular  and  irregular  waves.  The  surge 
notion  is  excluded  due  to  its  minor  practical  significance. 


*A  complete  listing  of  references  is  given  on  page  31. 
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Although  the  existing  analytical  method  has  its  own  merits  and  serves 

important  purposes,  it  has  been  felt  that  it  is  too  cumbersome  to  be  used 

in  the  conceptual  design  stage  at  which  numerous  candidate  hull  forms  are 
3 

examined.  Dalzell  has  addressed  this  problem  by  developing  a simplified 

method  of  computing  the  hydrodynamic  coefficients  associated  with  heave- 

pitch  coupled  motion  of  SWATH  ships  at  zero  speed  in  head  waves.  Since 

2 

the  major  portion  of  the  computer  program  for  prediction  of  motion  is  com- 
prised of  computing  the  sectional  added  mass  and  damping,  the  accomplish- 
ment of  Dalzell  in  reducing  the  size  and  computation  time  of  the  computer 
program  is  regarded  to  be  very  worthwhile. 

In  this  report,  an  extension  of  the  Dalzell  effort  to  all  hydrodynamic 
coefficients  associated  with  the  five  degrees-of-f reedom  of  motion  of  SWATH 
ships  is  described.  The  approach  employed  here  is  similar  to  the  Dalzell 
approach.  That  is,  the  sectional  wavemaking  damping  coefficients  of  SWATH 

cross  sections  are  obtained  using  the  far-field  potential-flow  theory.  By 

4 

application  of  the  Kramers-Kronig  relations  to  the  damping  coefficients, 
the  added-mass  coefficient,  which  can  be  regarded  as  the  conjugate  pair  of 
the  damping,  is  determined. 

The  advantage  of  this  procedure  over  the  conventional  multiple- 
expansion  method"*  or  source-distribution  method**  is  that  the  wavemaking 
damping  can  be  obtained  by  determining  the  outgoing  wave  amplitude  at  a 
far  field  rather  than  by  determining  the  pressure  distribution  on  the  body. 
The  outgoing  wave  amplitude  can  be  closely  approximated  by  representing  the 
flow  disturbances  by  a pulsating  single  source,  a dipole,  or  a combination 
of  both.  Furthermore,  the  Kramers-Kronig  relations  can  be  expressed  in 
terms  of  a half-range  double  Fourier  transform.  Because  an  efficient 
numerical  procedure  such  as  fast  Fourier  transform  (FFT)^  is  a commonly 
available  computer  routine,  the  process  of  obtaining  the  added-mass  coeffi- 
cients can  be  quickly  performed. 

Once  the  sectional  added  mass  and  damping  coefficients  in  heave,  sway, 
roll,  and  roll-sway  coupled  modes  are  obtained,  the  rest  of  the  motion  pre- 
diction process  follows  identically  to  that  described  in  Reference  1. 


2 


Comparison  of  the  results  obtained  by  experiments,  the  method  of 

g 

source  distribution,  and  the  present  approximate  method  are  presented  in 
graphs.  Agreements  are,  in  general,  fair.  The  approximate  method  yields 
'•  ’vings  in  the  computation  time  by  an  order  of  magnitude  compared  to  the 
conventional  method  of  Reference  8.  The  real  merit  of  the  approximate 
method  will  be  evaluated  in  the  future  when  an  analytical  method  of  pre- 
diction of  SWATH  motions  in  waves  is  developed. 


ANALYTICAL  METHODS 

BACKGROUND 

Within  a linear  analysis,  the  added  mass  and  damping  coefficients  of 
oscillating  twin-cylindrical  bodies  of  arbitrary  cross  section  have  been 

g 

determined  by  the  method  of  pulsating  source  distribution.  This  method 
has  shown  satisfactory  agreement  with  the  experimental  results  of  various 

g 

cross  sectional  shapes  such  as  twin  semicircles,  rectangles,  triangles, 

9 

and  SWATH  sections. 

To  obtain  reasonably  accurate  added  mass  and  damping  coefficients  of 
a SWATH  cross  section,  at  least  15  to  20  points  of  source  distribution  on 
the  submerged  contour  of  one  hull  should  be  taken  for  each  frequency  of 
oscillation.  To  obtain  all  the  hydrodynamic  coefficients  in  the  linear 
equations  of  motion  covering  all  degrees-of-f reedom,  the  sectional  added 
mass  and  damping  coefficients  should  be  computed  for  the  sway,  heave,  and 
roll  modes  for  about  20  cross  sections  along  the  length  of  a ship.  If  the 
motion  of  a SWATH  ship  in  irregular  waves  for  different  wave  headings  and 
ship  speeds  is  to  be  computed,  the  foregoing  number  of  calculations  should 
be  repeated  for  each  wave  heading,  ship  speed,  and  frequency  of  oscillation. 
Thus,  one  can  easily  expect  thac  a significant  amount  of  computer  time  will 
be  required  to  evaluate  the  seakeeping  qualities  of  a ship. 

In  spite  of  this  large  expense  of  computation  time,  what  we  are  ob- 
taining, at  best,  are  motions  based  on  the  hydrodynamic  coefficients  which 
are  obtained  under  the  assumption  of  two-dimensional  flow  conditions  at 
each  cross  section.  Furthermore,  unlike  monohull  ships,  the  motion  of 
SWATH  ships  cannot  be  predicted  accurately  by  using  only  the  coefficients 


obtained  from  the  potential-flow  theory.  There  are  viscous  effects  and 
fin  effects  to  be  included  in  the  hydrodynamic  coefficients.  It  seems 
reasonable,  therefore,  to  seek  a simplified  method  of  obtaining  the  sec- 
tional hydrodynamic  coefficients,  if  the  economic  savings  in  the  labor  and 
time  can  justify  the  loss  of  accuracy  of  the  results  within  a certain 
margin. 

ANALYSIS 

From  the  principle  of  conservation  of  energy,  the  wavemaking  damping 
of  an  oscillating  two-dimensional  body  can  be  determined  by  the  amplitude 
of  the  radiating  waves  at  a far  distance  from  the  body.  This  implies  that 
the  damping  coefficient  can  be  determined  by  examining  the  far-field  be- 
havior of  the  disturbances  generated  by  the  oscillating  body.  Once  the 
wavemaking  damping  is  known,  the  so-called  Kramers-Kronig  relations  can  be 
invoked  to  obtain  the  added  mass  of  the  body  Thus,  both  coefficients  can 
be  obtained  by  examining  the  far-field  behavior  of  the  outgoing  waves. 

The  advantage  of  examining  the  far-field  behavior  of  the  fluid  disturbances 
lies  in  the  fact  that  the  solution  of  the  potential  function  can  be  ob- 
tained by  examining  the  asymptotic  behavior  of  the  function  as  the  hori- 
zontal coordinate,  for  example  y,  goes  to  infinity. 

The  analysis  will  be  carried  out  in  the  following  manner.  First,  the 
potential  function  representing  the  outgoing  wave  at  the  far  field,  which 
is  generated  by  oscillating  the  demihull  of  a two-dimensional  SWATH  form, 
will  be  determined.  Then,  the  interference  and  blocking  effects  of  the 
other  hull  on  the  outgoing  waves  at  the  far  field  will  be  determined.  The 
next  step  is  to  determine  the  damping  coefficient  in  terms  of  the  wave 
amplitude  and  then,  the  added-mass  coefficient  by  the  Kramers-Kronig 
relations. 

Heave  Damping  Approximation 

Assume  that  the  circular  portion  of  a cross  section  of  the  demihull 
of  a SWATH  ship  can  be  represented  by  a vertical  dipole,  with  the  pulsating 


4 


velocity  V located  at  a certain  point  on  the  submerged  area  of  the  demihull 
cross  section,  which  is  shown  in  Figure  1.  The  notations  denoting  the 
dimensions  of  the  cross  section  are  also  given  in  Figure  1. 

The  velocity  potential  representing  the  vertical  dipole  in  an  un- 
bounded fluid  is  given  by 


♦d--is 


33  | cos  6 
P r 


where  p = water  density 

= submerged  cross  sectional  area  of  the  demihull 

= heave  added  mass  of  the  demihull  cross  section  in  an 
unbounded  fluid 

(r,0)  = cylindrical  coordinates  such  that  y = r sin  8 and 

z = r cos  0 - d , where  d is  the  depth  of  the  point  where 
o o 

the  dipole  is  located 


The  heave  added  mass  be  approximated  by 


W) 


Due  to  the  presence  of  the  surface-piercing  vertical  strut,  the  dis- 
placed area  of  the  cross  section  changes  according  to  the  vertical  oscilla- 
tion of  the  body.  The  change  in  the  displaced  area  is  equivalent  to  the 
flux  through  the  width  of  the  strut  which  is  Vb^.  Thus,  an  equivalent 
source  strength  to  generate  this  flux  can  be  obtained  by  Vbo/(2ir),  and  the 
corresponding  source  located  at  the  point  (0,-dg)  is  expressed  by 


2Tln  |/y2+(z+ds)2 


Figure  1 - Demihull  Cross  Section  of  a SWATH  Ship 


Thus,  the  velocity  potential  representing  a heaving  demihull  SWATH 
cross  section  is  obtained  approximately  by 


= ^d  + 


S 


v_ 

2tt 


a Vb 

cos  9 . o „ 

+ ~^r 

r 2tt 


2 2 
y +(z+d s) 


(4) 


The  far-field  behavior  of  the  singularities  represented  by  Equation 

(4),  where  the  z = 0 is  regarded  as  the  calm-water  line,  can  be  obtained 

, 11 
by 


-(v£) 


KeK(z-do+iy) 


+ b eK(z-ds+iy) 
o 


(-itoh  e ) 
o 


(5) 


where  only  the  real  part  of  the  right-hand  side  is  to  be  realized,  and 

„ . , -iwt 

V = -icoh  e 
o 

2 

K = 0)  /g,  the  wave  number 
0)  = circular  frequency  of  oscillation 
g = gravitational  acceleration 
i = 

Hereafter,  when  a real  function  is  expressed  in  a complex  form  involving 
e it  will  be  tacitly  assumed  that  only  the  real  part  of  the  complex 

expression  is  meant. 

Because  the  wave  amplitude,  which  is  denoted  by  A^,  is  obtained  by 


A 


1 


g z = 0 

I y I = 00 


then 
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4-k(v^) 


+b  Ke 
o 


From  the  conservation  of  energy,  the  rate  of  work  imparted  on  the 
fluid  by  the  body  during  one  cycle  of  oscillation  should  be  equal  to  the 
rate  of  energy  carried  by  the  radiating  surface  waves  during  the  same 
period.  This  means  that 


T 0 


J F.  V dt  = -2p  J dt  J <ft<p 


where 


= period  of  oscillation 


= f-u)  (M+a,,)  cos  u)t-u)b..  sin  wt+c . . cos  o)t]h 
jj  jj  JJ  ° 


-(oh  sin  wt 
o 


gA^  cos  (Ky-u)t+a)e 


= phase  angle  with  respect  to  the  vertical  motion  of 
the  body 


(a. . ,b. . ,c, ,)  = (added  mass,  damping,  restoring)  in  the  jth  mode 
^ due  to  the  motion  in  the  jth  mode 


= mass  of  the  body 


From  Equation  (7) 


bjj  ' p" 


(y 


Substituting  Equation  (6)  into  (8) 


b3‘x>  - p» 


where  the  superscript  1 is  used  to  indicate  the  quantities  corresponding 
to  the  demihull  only. 
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To  take  into  account  the  interference  or  blockage  effects  from  the 
other  hull,  the  reflected  images  at  y = -2b  of  the  dipole  and  source  in- 
troduced in  the  foregoing  analysis  should  be  considered.  Further,  assume 
that  the  demihull  is  replaced  by  a vertical  barrier  having  the  depth  D; 
hence,  the  waves  generated  by  one  hull  cannot  be  transmitted  completely 
beyond  the  other  hull. 

If  we  let  £ denote  the  ratio  of  the  wave  amplitude  of  the  incident 
wave  before  the  barrier  to  the  amplitude  of  the  transmitted  wave  behind 
the  barrier,  the  waves  generated  by  the  singularities  can  be  expressed  in 
the  form 

£ = cos  (Ky-wt)  + E cos  (Ky+2b-<jJt+aT)  (11) 

where  the  incident  wave  amplitude  is  assumed  to  be  unity,  and  a^,  is  the 
change  of  the  phase  angle  of  the  transmitted  wave. 

Equation  (11)  can  be  expressed  as 

£ = Bq  cos  (Ky-u)t+6)  (12) 

where 


(13) 


Bq  = [1+e2+2e  cos  (2Kb+aT)  ]1//2 


6 = tan 


^ -e  sin  (2Kb+aT) 
1+E  cos  (2Kb+a^,) 


The  transmission  coefficient  E and  the  phase  angle  oi  are  given  by^ 


Kx(KD) 

/,tt2I12(KD)+K12(KD) 


aT  = tan 


_1/Tai(KD)\ 
\ k2 (KD)  ) 


where  1^  and  are  the  modified  Bessel  functions.  For  x < 1,  it  can  be 
shown  that 

..  3 5 7 


[ (v)  = X + X.  + X + 

1 2 16  384  18432 


KX(X) 


x + ln(i)I1<x) 


(#»  3 (?-») 


(1  -»)  , 

36864  X ‘ 


where 


Y = Him 


= 0.57721566. 


is  known  as  the  Euler  constant.  For  KD  _>  1,  it  can  be  approximated  that 
e = 0 and  aT  = 0. 
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Because  3q  represents  the  interference  and  blockage  effects  on  the 
wave  amplitude  at  the  far  field,  the  twin-hull  damping  coefficient  can  be 
determined  by 


b 


33 


r E (Dr  2 
C1  b33  6o 


(18) 


where  is  introduced  here  as  a correction  factor  if  the  need  arises. 

When  the  strut  thickness  is  zero,  the  source  potential  given  by 
Equation  (3)  can  be  disregarded  because  bQ  = 0. 

Also,  due  to  the  absence  of  the  vertical  barrier,  it  is  assumed  that 
e = 1 and  = 0.  Thus,  the  heave  damping  of  completely  submerged  twin 
circles  can  be  obtained  from  Equations  (9)  and  (18)  as 
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where 


S = ttR2 
A 


dc  = the  depth  to  the  center  of  the  circle 


(19) 
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The  last  two  terms  in  the  bracket  of  the  expression  of  11^3  represent,  re- 
spectively, the  twin-hull  effect  and  the  free-surface  effect  which  is  as- 
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sumed  to  be  a rigid  wall. 

Sway  Damping  Approximation 

From  a far-field  view,  the  demihull  of  a SWATH  section  in  the  sway 
motion  can  be  regarded  as  a vertical  plate  of  the  same  draft  as  the  body 
draft.  For  small  frequencies  of  oscillation  which  are  of  practical 
interest  to  SWATH  motion,  it  can  be  further  assumed  that,  near  the  body, 
the  fluid  disturbances  would  be  close  to  those  made  by  a vertical  plate  of 
twice  the  length  of  the  body  draft  in  sway  motion  in  an  unbounded  fluid. 

If  the  sway  velocity  of  the  body  motion  is  denoted  by  U,  the  flow 
field  can  be  represented  by  the  velocity  potential  which  can  be  expressed 
in  terms  of  the  dipole  distribution  as 


(20) 


1/2  2 

where  the  dipole  density  \J  D -£  can  be  obtained  from  the  well-known  solu- 
tion of  a flat  plate  subject  to  normal  uniform  flow  by  taking  the  potential 
jump  across  the  plate. 

The  corresponding  far-field  expression  which  represents  the  radiating 
free-surface  waves  is  obtained  from  Equation  (13.29)  of  Reference  11  as 
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where 
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in  which  is  the  modified  Struve  function. 


The  amplitude  of  radiating  wave  at  |y|  = 00  is  obtained  by 


2o)  1 U | K£  (K) 
g 


and  if 


then 


U = -iu)h  e 
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-iojt 


A 2K2B(K) 


(23) 


Thus,  the  sway  damping  for  demihull  of  a SWATH  cross  section  can  be  ob- 
tained from  Equations  (8)  and  (23)  as 


s a>  = -X- 
22  P“SA  l K2  ) 


4K2B2(K) 

SA 


(24) 


The  interference  effect  from  the  other  hull  will  be  treated  the  same 
as  in  the  previous  case.  That  is,  the  far-field  wave  amplitude  A^  will  be 
modified  as  3qA^  where  8Q  is  given  by  Equation  (13).  Consequently,  the 
sway-damping  coefficient  for  the  twin  hulls  is  obtained  by 
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where  C2  is  a correction  factor  which  should  be  determined  as  the  need 
arises. 

For  the  cross  section  without  the  struts,  the  sway  damping  is  assumed 
to  be  similar  to  the  heave  damping  given  by  Equation  (19),  i.e., 
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M22 

— — - = one-half  of  the  sway  added  mass  of  twin  circles  under  a 


rigid  surface 

- 
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Roll  Damping  Approximation 

The  roll  moment  exerted  on  a SWATH  cross  section  due  to  hydrodynamic 
pressures  p can  be  expressed  by 


^ = j)  p(ynz-zny)d£ 


where  n and  n are,  respectively,  the  y-component  and  the  z-component  of 

y z 

the  unit  normal  vector  pointing  into  the  body  on  the  submerged  contour  of 
the  cross  section,  and  £ di  is  the  integral  over  the  submerged  cross- 
section  contour. 
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If  <|>R  denotes  the  velocity  potential  associated  with  a forced  roll 
oscillation  of  the  cross  section,  then,  from  the  linearized  Bernoulli 
equation. 


P - ipoxj) 


hence 


Mj^  = ipoj  ^ <f>R(ynz-zn  )d«. 


Assume  that 
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Then,  substituting  Equation  (29)  into  Equation  (28), 
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where  the  mean-value  theorem  is  used  by  defining 
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The  last  integral  in  Equation  (30)  can  be  approximated  by 


2 2 2 

where  it  is  assumed  that  y = b . An  approximate  z is  given  by 
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For  the  cross  section  without  the  struts,  the  roll  damping  is  approxi- 
mated by 


b44  = b33  l1+  7T 

D 


(33) 


where  b^  corresponds  to  the  one  given  by  Equation  (19). 


Sway-Roll  Coupling  Damping  Approximation 

The  same  assumptions  made  in  the  approximation  of  the  roll  damping 
will  yield 
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For  the  cross  section  without  the  struts, 
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where  b ^ is  given  by  Equation  (26). 
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Added  Mass  Approximation 

Once  the  damping  coefficients  are  determined,  then  the  corresponding 
conjugate  pairs,  i.e.,  the  added-mass  coefficients  can  be  obtained  by  the 
Kramers-Kronig  relations.  A convenient  form  for  numerical  manipulations 
is  given  in  Reference  4 as 


a (co)  - a..  (°°)  = 
jk  ' 


cos  cot 


r bik((o,)"bik(oo) 

J Sin  U’t  dco'dt 


for  j , k = 2,  3,  and  4 

Wh6re  3jk  ls  the  added-mass  coefficient  in  the  jth  mode  due  to  the  motion 
in  the  kth  mode,  and  bjk(«),  in  general,  is  zero  because  no  surface  wave 
can  be  generated  at  w = The  nondimens ional  added-mass  coefficients  will 
be  denoted  by  the  bar  sign,  and  are  defined  by 

(a  5 j = (a22>a33} 

U22’a33'  ~ pS  (38) 
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£ = £ - S24 
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The  added-mass  coefficients  at  the  infinite  frequency  are  approximated 


i (oo)  = 1D_ 
a22 ' ’ ttS, 


which  is  based  on  the  flat  plate  results, 
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where  Tr  D / 32  is  obtained  by  a,  , (°°) /a„? (°°)  for  a vertical  plate  having  the 
14  44 

draft  D, 
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For  the  cross  section  without  the  struts,  assume  that 
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a42 (00)  = a22  r 
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where  is  a correction  factor  which  should  be  determined  when  the  need 
arises.  The  expressions  given  by  Equations  (45)  and  (46)  correspond  to 
the  cases  that  twin  circles,  separated  by  2b  between  the  centers,  move  with 
a constant  velocity  in  the  unbounded  fluid  in  the  direction  opposite  to 
each  other,  and  in  the  same  direction  normal  to  the  line  connecting  the 
centers  of  the  circles,  respectively. 
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DISCUSSION  OF  RESULTS 


The  numerical  results  obtained  by  the  approximate  methods  described  in 
the  preceding  sections  are  presented  in  Figures  2 through  9.  The  results 
of  the  approximate  methods  are  compared  with  the  results  obtained  by  the 
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accurate  theory,  and  with  the  available  experimental  results.  Hereafter, 
the  method  described  in  Reference  8 shall  be  referred  to  as  "accurate 
theory. " 

The  damping  coefficients  of  a heaving,  two-dimensional  SWATH  cross 
section  are  presented  in  Figure  2.  The  heave  damping  b^  is  made  non- 
dimensional  by  the  product  of  the  displaced  fluid  mass  (pS^)  and  the  cir- 
cular frequency  of  oscillation  (ui) . The  relative  dimensions  of  the  section 
are  as  shown  in  the  figure.  The  approximate  heave  damping  coefficients  are 
obtained  by  Equation  (18)  with  = 1 and  aT  = 0.  It  means  that  the  phase 
change  of  the  transmitted  wave  beyond  the  vertical  barrier  is  neglected. 
Inclusion  of  the  phase  change  in  Equation  (18)  resulted  in  poorer  agree- 
ment with  the  accurate  theoretical  results.  As  can  be  seen,  the  three 
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results  agree  reasonably  well  for  the  frequency  number  (uj  R/g)  less  than 
0.5. 

Because  the  shortest  wavelength  of  practical  interest  in  SWATH  motion 
at  zero  forward  speed  is  about  one-half  of  the  ship  length,  the  highest 
frequency  number  corresponding  to  this  wavelength  for  a SWATH  ship  having 
the  length-to-diameter  ratio  of  15  is  about  0.8.  Beyond  this  highest 
frequency  limit,  SWATH  ships  would  barely  respond  to  waves,  which  means 
that  there  is  no  need  of  computing  the  transfer  function  of  the  motion 
response  beyond  this  frequency  limit.  Thus,  the  approximate  results  shown 
in  Figure  2 would  be  quite  reasonable  for  use  in  the  computation  of  motion 
of  SWATH  ships  in  waves. 

The  heave  added-mass  coefficient  for  the  SWATH  section  shown  in 
Figure  2 is  presented  in  Figure  3.  The  results  from  the  approximate  method 
are  obtained  from  Equations  (37)  and  (42).  The  integrals  on  the  right-hand 
side  of  Equation  (37),  in  fact,  represent  double  Fourier  transforms.  A 
fast  Fourier  transform  subroutine  of  the  system  library  of  the  computer 
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Figure  3 - Heave  Added-Mass  Coefficient  of  a SWATH  Section 


facility  at  the  Center  was  used  to  evaluate  the  double  Fourier  transforms. 
The  subroutine  is  called  "FFT"  which  was  originally  developed  at  the  Los 
Alamos  Scientific  Laboratory  in  California. 

As  can  be  observed  in  Figure  3,  an  apparent  discrepancy  between  the 
exact  and  approximate  methods  is  at  those  frequencies  at  which  the  exact 
method  shows  abrupt  discontinuities.  So  far  as  the  problem  is  confined  to 
,1  two-dimensional  flow  condition,  the  discontinuities  exhibited  by  the  exact 
method  can  be  real  physical  phenomena*;  although,  the  experimental  results, 
unfortunately,  were  not  taken  at  those  frequencies.  However,  for  three- 
dimensional  twin  bodies  such  as  SWATH,  there  cannot  be  a complete  trapping 
of  waves  between  the  two  hulls;  hence,  no  abrupt  discontinuities  of  the 
hydrodynamic  coefficients  at  certain  frequencies  should  exist.  In  this 
respect,  the  approximate  method,  which  does  not  create  the  discontinuities, 
may  cause  less  errors  at  those  frequencies  at  which  the  discontinuity  in 
the  hydrodynamic  coefficient  occurs. 

The  SWATH  cross  sections  at  far  forward  or  aft  portion  of  the  body  do 

not  have  vertical  struts.  For  these  sections,  the  cross-section  view  is 

completely  submerged  twin  circles.  The  heave  added  mass  and  and  damping 

coefficients  of  such  a cross  section  are  shown  in  Figure  4.  For  the 
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damping,  Equation  (19)  was  used  with  = e , and  for  the  added  mass. 
Equations  (37)  and  (46)  with  = 1.1  were  used.  Agreement  between  the 
two  methods  appears  good  from  a practical  viewpoint  for  the  prediction  of 
SWATH  motion.  This  is  because  a significantly  larger  magnitude  of  heave 
damping  contributed  by  viscosity  and  stabilizing  fins  should  be  added  to 
the  wavemaking  damping  in  the  analytical  prediction  of  SWATH  motion;  there- 
fore, a 30  p-jrcent  error,  as  can  be  observed  in  the  heave  damping  in  Figure 
4,  is  not  going  to  affect,  significantly,  the  final  results  of  the  motion. 
There  appears  to  be  a maximum  of  about  5 percent  error  in  the  heave  added 
mass  calculated  by  the  approximate  method.  Because  the  SWATH  motion  is  far 
less  sensitive  to  the  added  mass  than  to  the  damping,  a 5 percent  error  in 
the  added-mass  coefficient  is  not  going  to  affect,  significantly,  the  pre- 
dicted results  of  SWATH  motion. 


*Such  evidence  is  clearly  demonstrated  by  the  experiments  for  several 
two-dimensional  twin  bodies,  as  described  in  Reference  8. 
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The  sway  damping  of  a SWATH  cross  section  is  presented  in  Figure  5, 
together  with  the  relative  dimension  of  the  cross  section.  The  damping  co- 
efficients obtained  by  the  approximate  method  are  based  on  Equation  (25) 

(dQ/(2R)-KD) 

with  the  modification  factor  C.  = e wher  d , the  depth  of  the 
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centroid  of  the  submerged  cross  sectionai  area,  can  be  obtained  by 

do  ■ ^ [l  bodS2+’R2ws+R)] 


At  present,  there  is  no  rational  explanation  for  the  derivation  of  the 
modification  factor  Cj  as  it  is  except  that  such  a factor  seems  to  repre- 
sent the  behavior  of  the  sway  damping  in  the  low  and  high  frequency  ranges 
reasonably  well.  Also,  in  the  vicinity  of  the  frequency  where  an  abrupt 
discontinuity  occurs,  the  modification  factor  appears  to  make  the  approxi- 
mate method  behave  smoothly  across  the  discontinuity  predicted  by  the 
accurate  theory.  The  chained  curve  shown  in  Figure  5 represents  the 
results  of  the  approximate  method  using  = 1. 

The  sway  added  mass  of  the  SWATH  cross  section  shown  in  Figure  5 is 
presented  in  Figure  6.  The  approximate  results  were  obtained  by  Equations 
(37)  and  (41).  In  the  frequency  range  of  0.2  < oo^R/g  < 0.5,  the  approxi- 
mate results  show  a large  discrepancy.  It  is  not  obvious  at  the  present 
stage  whether  the  sway  added  mass  computed  by  the  approximate  method  would 
necessarily  result  in  poorer  prediction  of  SWATH  motion  than  that  obtained 
by  the  accurate  theory. 

The  sway  added  mass  and  damping  coefficients  for  submerged  twin  circles 
are  presented  in  Figure  7.  Equation  (26)  was  used  to  obtain  the  damping, 
and  Equations  (37)  and  (45)  with  = 1.1  were  used  to  obtain  the  added 
mass.  The  trend  of  the  results  is  very  close  to  that  of  the  heave  added 
mass  shown  in  Figure  4. 

The  roll  damping  of  a SWATH  cross  section  is  presented  in  Figure  8. 

The  dimensions  of  the  cross  section  are  also  shown  in  the  figure.  Note 
that  the  center  of  roll  is  located  above  the  calm  water  level  by  the 
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Figure  7 - Sway  Added-Mass  and  Damping  Coefficients  of 
Submerged  Twin  Circles 


distance  of  3.1  times  the  radius  of  the  cross  section.  The  damping  coeffi- 
cients were  obtained  by  Equation  (31).  Both  the  experimental  and  the 
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accurate  theoretical  results  show  discontinuity  at  about  u)  R/g  = 0.65. 
However,  as  in  the  case  of  heave  and  sway,  the  approximate  results  do  not 
exhibit  any  discontinuity.  The  relative  trends  of  the  three  different 
results  are  quite  similar  to  the  sway  damping  shown  in  Figure  5. 

The  roll  added  inertia  coefficient  of  the  cross  section  shown  in 

Figure  8 is  presented  in  Figure  9.  The  approximate  results  show  a large 
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discrepancy  in  0.2  < a)  R/g  < 0.5  as  was  already  observed  in  Figure  6 for 
the  sway  added  mass.  The  same  remarks  as  given  earlier  on  the  sway  added 
mass  apply  to  the  roll  added  inertia  obtained  by  the  approximate  method. 

The  accurate  theory,  referred  to  frequently  in  this  report,  is  based 
on  the  Green  function  method  which  leads  to  an  integral  equation  involving 
the  unknown  strengths  of  the  source  distribution  over  the  submerged  contours 
of  the  cross  section  of  the  twin  hulls.  The  numerical  evaluation  of  the 
integral  equation  requires  a segmentation  of  the  integral  along  the  cross- 
section  contour  into  a finite  number  of  integrals.  Increasing  the  number 
of  segments  of  the  contour  increases  the  accuracy  of  the  results.  For 
SWATH  sections,  it  takes  about  fifteen  to  twenty  segments  on  the  contour  of 
the  demihull  cross  section  to  obtain  accurate  results.  A typical  computa- 
tion by  the  CDC  6000-series  computer  at  the  Center  for  the  sway  and  heave 
hydrodynamic  coefficients  for  submerged  twin  circles  for  twenty  frequencies 
took  about  220  seconds  of  executing  time. 

On  the  other  hand,  the  present  approximate  method  took  only  about  16 
seconds  to  obtain  the  sway,  heave,  roll,  and  sway-roll  coupling  coefficients 
for  256  frequencies.*  Even  if  a discount  is  made  for  the  large  number  of 
frequencies  which  are  dictated  by  the  usage  of  the  fast  Fourier  transform 
technique  but  not  by  the  practical  necessity,  the  approximate  theory  seems 
to  yield  a time  saving  in  the  computation  by  an  order  of  magnitude  compared 
to  the  accurate  theory. 


*Most  of  the  fast  Fourier  transform  subroutine  requires  that  the  func- 
tion to  be  transformed  should  be  given  in  numbers  of  2n  where  n is  an 
integer.  For  the  present  work,  it  was  found  that  n = 8 yields  satisfactory 
results. 
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Figure  9 - Roll  Added-Inertia  Coefficient  of  a SWATH  Section 
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SUMMARY  AND  CONCLUSIONS 

To  develop  a simplified  analytical  prediction  of  SWATH  motions  in 
waves,  approximate  methods  of  computing  the  hydrodynamic  coefficients  of 
oscillating  cross  sections  of  a SWATH  configuration  are  developed.  The 
hydrodynamic  coefficients  studied  are  the  added  mass  and  the  wavemaking 
damping  coefficients  in  sway,  heave,  roll,  and  sway-roll  coupling  modes. 

The  approximate  methods  are  based  on  the  concept  that  the  wavemaking 
damping  coefficients  can  be  obtained  by  the  amplitudes  of  the  outgoing 
waves  at  the  far  field.  The  outgoing  waves  at  the  far  field  can  be  repre- 
sented by  a pair  of  sources  and  of  dipoles  for  the  heave  oscillation  and 
by  a dipole  distribution  on  the  vertical  centerline  of  each  hull  for  the 
sway  and  roll  oscillations.  The  strength  of  the  singularities  are  obtained 
under  the  assumption  that  the  cross  section  is  oscillating  in  an  unbounded 
fluid. 

The  hydrodynamic  interactions  between  the  two  hulls  are  treated  as  if 
one  hull  is  replaced  by  a vertical  barrier  having  the  draft  as  the 
original  form  and,  therefore,  it  blocks  the  passage  of  the  waves  generated 
by  the  other  hull. 

The  added-mass  coefficients  are  obtained  by  invoking  the  Kramer s- 
Kronig  relations  with  the  known  damping  coefficients  and  the  added-mass 
coefficients  at  the  infinite  frequency,  as  indicated  by  Equation  (37).  The 
evaluation  of  Equation  (37)  is  carried  out  by  a well-known  numerical 
algorithm  called  fast  Fourier  transform  technique.  When  needed,  the  ap- 
proximate methods  are  augmented  by  the  modification  factors  which  are  de- 
termined by  a numerical  trial-and-error  approach. 

The  results  from  the  approximate  methods,  the  accurate  theory  based 
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on  the  source  distribution  on  the  submerged  contours,  and  the  experiments, 
if  available,  are  compared  in  Figures  2 to  9.  The  agreements  of  the 
results  in  the  frequency  range  of  practical  interest  for  SWATH  motion  are, 
in  general,  fair;  however,  the  abrupt  discontinuities  of  the  hydrodynamic 
coefficients  at  certain  frequencies  exhibited  by  the  accurate  theory  and 
the  experiments  are  not  reproduced  by  the  approximate  methods.  The  abrupt 
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discontinuities  are  merely  caused  by  the  two-dimensional  flow  condition, 
and  in  reality,  for  three-dimensional  bodies,  they  do  not  exist;  hence, 
the  continuity  of  the  hydrodynamic  coefficients  at  all  frequencies  given  by 
the  approximate  methods  could  be  more  realistic. 

The  time  saving  in  the  computation  of  the  hydrodynamic  coefficients 
achieved  by  use  of  the  approximate  method  is  remarkable;  the  approximate 
method  can  provide  the  results  in  less  than  one-hundredth  of  the  time  which 
would  be  required  by  the  accurate  theory. 

It  is  too  early  to  judge  whether  the  approximate  methods  developed 
here  will  provide  reasonable  prediction  of  SWATH  motions  in  waves.  The 
usefulness  of  the  approximate  methods  will  be  realized  only  when  they  are 
shown  to  provide  adequate  predictions  of  the  motions  for  use  during  the 
stage  of  conceptual  hull  design  of  SWATH  configurations.  However,  it  is 
believed  that  a final  useful  method  of  predicting  SWATH  motion  can  be  de- 
veloped on  the  basis  of  the  approximate  methods  presented  in  this  report  by 
an  intelligent  usage  of  the  modification  factors  associated  with  the  ap- 
proximate methods. 
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